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Abstract 

Smart  material  transducers  employing  piezoceramic 
or  magnetostrictive  drive  components  typically  exhibit 
constitutive  nonlinearities  and  hysteresis  at  moderate  to 
high  drive  levels.  While  feedback  mechanisms  or  careful 
choice  of  operating  regimes  can  often  reduce  these  ef¬ 
fects,  spillover  into  high  frequency  dynamics  and  phase 
lags  associated  with  the  two  phenomena  will  degrade 
the  controller  performance  at  high  drive  levels  in  the 
absence  of  additional  compensation  or  nonlinear  con¬ 
trol  design.  In  this  paper,  we  discuss  two  techniques 
to  compensate  for  hysteresis  in  high  performance  trans¬ 
ducers.  The  first  is  based  on  a  complete  transducer 
model,  and  the  resulting  compensator  accommodates 
both  the  constitutive  nonlinearities  and  hysteresis  in¬ 
herent  to  the  smart  material  components.  The  second 
technique  employs  a  partial  inverse  compensator  based 
on  anhysteretic  models  for  the  material  behavior.  This 
accommodates  the  constitutive  nonlinearities  but  does 
not  incorporate  the  hysteresis;  the  latter  phenomenon  is 
then  addressed  through  the  inclusion  of  a  feedback  loop 
in  the  controller.  The  performance  of  the  partial  inverse 
compensator  is  illustrated  in  the  context  of  a  high  force 
Terfenol-D  transducer. 

1.  Introduction 

Smart  material  transducers  employing  piezoceramic, 
electrostrictive  or  magnetostrictive  cores  provide  the  ca¬ 
pability  for  both  actuating  and  sensing  in  a  wide  range  of 
high  performance  applications.  Magnetostrictive  trans¬ 
ducers  provide  large  forces  with  broader  bandwidths 
than  traditional  magnetic  shakers.  Piezoceramic  and 
electrostrictive  patches  have  the  capability  for  sensing 


and  actuating  while  minimally  affecting  the  passive  dy¬ 
namics  of  the  underlying  structure.  Furthermore,  they 
can  be  configured  in  a  variety  of  geometries  and  are  rela¬ 
tively  inexpensive  to  manufacture.  In  combination  with 
their  dual  sensing  and  actuating  capabilities,  these  prop¬ 
erties  make  them  attractive  in  a  number  of  control  sys¬ 
tems.  In  their  biased  state,  all  of  these  materials  exhibit 
an  approximately  linear  behavior  at  low  drive  levels. 
At  moderate  to  high  drive  levels,  however,  all  exhibit 
constitutive  nonlinearities  including  saturation  effects. 
Furthermore,  piezoceramics  and  magnetostrictives  ex¬ 
hibit  significant  hysteresis  while  relaxor  ferroelectrics  are 
hysteretic  when  employed  in  their  ferroelectric  phase. 
These  constitutive  nonlinearities  and  hysteresis  must  be 
accommodated  to  achieve  the  full  capabilities  of  the  ma¬ 
terials  in  high  performance  control  design. 

A  common  strategy  is  to  limit  input  levels  to  ranges 
in  which  the  material  behavior  is  approximately  linear. 
This  permits  the  use  of  linear  models  and  control  de¬ 
signs.  While  feasible  in  certain  applications,  this  sig¬ 
nificantly  limits  the  capability  of  the  transducers.  For 
piezoceramic  materials,  current  control  in  lieu  of  volt¬ 
age  control  can  be  used  to  maintain  an  approximately 
linear  response  [3].  However,  these  techniques  do  not 
appear  directly  applicable  to  magnetostrictive  transduc¬ 
ers.  For  certain  applications,  feedback  loops  can  be  used 
to  reduce  the  effects  of  hysteresis.  For  high  accuracy  or 
micropositioning  applications,  however,  further  compen¬ 
sation  is  typically  required.  This  can  be  achieved  either 
through  the  development  of  inverse  compensators  incor¬ 
porated  in  linear  control  designs  or  nonlinear  controllers. 
In  this  paper,  we  focus  on  the  development  of  full  or  par¬ 
tial  inverse  compensators  and  refer  the  reader  to  [4]  for 
an  example  of  nonlinear  optimal  control  design  in  the 
context  of  a  Terfenol-D  magnetostrictive  transducer. 
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The  concepts  underlying  full  and  partial  inverse  com¬ 
pensation  are  illustrated  in  Figure  1.  In  both  cases, 
we  assume  that  the  smart  material  transducer  provides 
hysteretic  inputs  to  the  system.  For  full  inverse  com¬ 
pensation,  a  complete  inverse  TV-1  is  constructed  to  ac¬ 
commodate  the  actuator  hysteresis  N  so  that  the  control 
input  u  to  the  plant  is  equal  (or  approximately  equal)  to 
the  desired  input  Ud  specified  by  the  control  law.  This 
inverse  filter  is  constructed  using  the  complete  model 
used  to  characterize  the  hysteretic  transducer.  The  par¬ 
tial  inverse  compensator  considered  here  accommodates 
the  constitutive  nonlinearities  through  the  anhysteretic 
curve.  The  uncompensated  hysteresis  then  causes  a 
phase  delay  between  u  and  ua  which  is  accommodated 
through  a  feedback  loop.  This  places  additional  impe¬ 
tus  on  the  control  design  but  significantly  reduces  the 
overhead  required  for  implementation. 

The  domain  wall  model  used  to  characterize  the  hys¬ 
teresis  is  summarized  in  Section  2.  This  model  is  ad¬ 
vantageous  in  that  it  is  equally  applicable  to  ferromag¬ 
netic  and  ferroelectric  materials  and  it  can  be  inverted 
to  provide  a  full  compensator  for  linear  control  design. 
Furthermore,  the  anhysteretic  component  provides  the 
filter  employed  in  partial  compensator  design.  The  de¬ 
velopment  of  partial  and  full  inverse  compensators  is 
summarized  in  Section  3.  The  performance  of  the  par¬ 
tial  inverse  compensator  is  illustrated  in  Section  4  for  a 
high  speed,  high  accuracy  milling  application. 

2.  Hysteresis  Model 


ferromagnetic  materials  throughout  the  remaining  dis¬ 
cussion.  The  reader  is  referred  to  [6]  for  an  analogous 
development  for  general  ferroelectric  materials  and  [7] 
for  details  regarding  the  use  of  the  model  to  character¬ 
ize  hysteresis  in  piezoceramic  materials. 

The  model  development  is  considered  in  two  steps.  In 
the  first,  the  anhysteretic  magnetization  Man  is  quan¬ 
tified  using  two  techniques;  the  first  is  based  on  Boltz¬ 
mann  statistics  while  a  polynomial  fit  comprises  the  sec¬ 
ond.  The  second  step  in  the  model  development  is  the 
quantification  of  the  irreversible  magnetization  Mjrr  and 
reversible  magnetization  Mrev  which  result  when  do¬ 
main  walls  respectively  translate  and  bend.  The  sum 
of  the  two  components  yields  the  total  magnetization. 

2.1.  Anhysteretic  Models 

The  anhysteretic  magnetization  at  a  point  in  the  ma¬ 
terial  is  dependent  on  the  effective  magnetic  field,  the 
saturation  magnetization  and  the  thermal  energy  of  the 
sample.  Under  the  assumption  of  constant  stress  ct0, 
the  effective  field  can  be  expressed  as  He  =  H  +  aM 
where  H  denotes  the  external  field  input  to  the  trans¬ 
ducer,  M  denotes  the  magnetization  and  the  parameter 
a  quantifies  the  effects  of  interdomain  coupling  and  mag¬ 
netoelastic  domain  interactions.  The  balance  of  thermal 
and  magnetostatic  energy  using  Boltzmann  principles 
respectively  yields  the  Langevin  and  Ising  spin  models 

Man  =  Ms£(He/a) 

Man  =  Ms  tanh  (He/a) 


The  model  described  here  is  based  on  the  tenet  that 
hysteresis  in  ferromagnetic  and  ferroelectric  materials  is 
due  to  the  impeded  movement  of  domain  walls  which 
are  pinned  at  inclusions  or  defects  in  the  material.  For 
brevity,  we  focus  on  on  the  modeling  of  hysteresis  in 


Figure  1.  (a)  Full  inverse  compensator  employed  as  a 
filter  for  linear  control  design;  (b)  Partial  inverse  com¬ 
pensator  based  on  the  anhysteretic  curve. 


for  the  anhysteretic  magnetization.  Here  the  Langevin 
function  is  defined  by  £{z)  =  coth(^)  —  1  /z.  The  pa¬ 
rameters  a  and  a  are  estimated  for  a  given  transducer 
through  a  least  squares  fit  to  experimental  data. 

A  second  technique  for  modeling  the  anhysteretic 
magnetization  is  through  the  use  of  piecewise  polyno¬ 
mials  or  general  power  law  representations  fit  to  data 
measured  in  specific  applications.  For  example,  a  piece- 
wise  fit  employing  three  components  has  the  form 


h(H)  , 

Hmin  <  H  <  H\ 

h{H)  , 

HxKH  <  h2 

(2) 

MH)  , 

H2  <  H  <  Hmax 

where  fi(H),  fe{H)  are  piecewise  functions  and 

Hmin,  Hmax  are  the  minimum  and  maximum  field  in¬ 
puts.  The  determination  of  specific  functions  to  attain 
an  appropriate  inverse  is  illustrated  in  Section  4. 

2.2.  Full  Hysteresis  Model 

As  detailed  in  [1],  the  irreversible  magnetization  is 
quantified  by  determining  the  energy  required  to  re¬ 
orient  dipoles  in  the  presence  of  an  applied  field.  This 
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yields  the  differential  equation 


d  j'  j' 

dH 


=  S 


Man 


kd  -  a  ( Man  -  Mirr) 


(3) 


where  the  parameter  k  provides  an  average  measure 
of  the  energy  required  to  re-orient  domains  and  hence 
translate  domain  walls.  The  parameter  6  is  defined  to 
have  the  value  +1  when  dH  >  0  and  —1  when  dH  <  0  to 
guarantee  that  pinning  always  opposes  changes  in  mag¬ 
netization.  In  applications,  d  can  be  directly  determined 
from  the  magnetic  field  data.  Finally,  the  parameter  <5 
is  defined  to  be  0  if  {dH  >  0  and  M  >  Man}  or  {dH  < 
0  and  M  <  Man }  and  is  1  otherwise  to  incorporate  the 
physical  observation  that  following  field  reversal  at  the 
loop  tip,  the  changes  in  magnetization  are  reversible  un¬ 
til  the  anhysteretic  value  is  reached  [2]. 

The  second  component  of  the  magnetization  is  the  re¬ 
versible  magnetization  which  reflects  the  degree  to  which 
domain  walls  bend  before  attaining  the  energy  necessary 
for  translation.  As  derived  in  [1],  the  inclusion  of  this 
component  yields  the  total  magnetization 


M  =  (1  —  c)Mirr  +  cMan  (4) 


where  Mjrr  is  defined  by  (3)  and  the  anhysteretic  mag¬ 
netization  Man  is  given  by  (1)  or  (2). 

To  provide  a  relation  which  facilitates  inversion,  it 
is  advantageous  to  express  the  output  magnetization  as 
a  function  of  the  input  field.  As  detailed  in  [5],  the 
magnetization  at  a  given  field  level  can  be  specified  as 
the  solution  to  the  differential  equation 


Sr  =  *"■") 


where  T  is  given  by 
T{H,  M)  = 


l+cMsa§§(«±^) 


(5) 


S  [MSL  _  M] 

kd  -q[Ms£(^M)  _  M] 


—cMs 


dC  ( H  +  aM ' 
dH 


2.3.  Full  Model  Inverse 

The  hysteresis  model  in  the  form  (5)  is  amenable  to 
inversion  through  consideration  of  the  complementary 
differential  equation 


dM-1  _  1 

dH  ~  T{M-',H) 


(7) 


with  T  defined  by  (6).  This  provides  the  exact  inverse 
N^1  if  the  parameters  Ms,a,  a,  c  and  k  are  known  ex¬ 
actly  or  an  approximate  inverse  IV-1  if  the  parameters 


are  unknown  or  slowly  changing  due  to  changing  oper¬ 
ating  conditions.  Details  regarding  the  design  of  control 
methods  which  employ  an  approximate  inverse  can  be 
found  in  [9]. 

3.  Full  and  Partial  Inverse  Compensators 

The  previously  described  anhysteretic  and  hysteresis 
models  provide  the  capability  for  constructing  both  full 
and  partial  inverse  compensators  for  linear  control  de¬ 
sign.  To  construct  a  full  inverse  compensator,  the  dif¬ 
ferential  equation  (7)  is  used  to  construct  a  filter  TV-1 
which  is  employed  in  the  manner  illustrated  in  Fig¬ 
ure  la.  This  provides  a  method  for  linearizing  the  sys¬ 
tem  when  the  degree  of  hysteresis  is  large.  Details  re¬ 
garding  the  performance  of  the  full  inverse  compensator 
for  a  Terfenol-D  magnetostrictive  transducer  are  pro¬ 
vided  in  [8].  We  note  that  for  the  full  inverse  to  be 
successful,  it  is  important  to  have  accurate  initial  con¬ 
ditions  which  is  difficult  to  attain  experimentally.  This 
is  one  motivation  for  considering  a  partial  compensator. 

A  partial  inverse  compensator  can  be  constructed  us¬ 
ing  the  Ising  spin  model  (1)  or  piecewise  power  repre¬ 
sentation  (2).  The  resulting  inverse  incorporates  the 
nonlinear  behavior,  including  saturation  effects,  but  ne¬ 
glects  the  phase  delays  due  to  hysteresis  (see  Figure  lb) . 
For  systems  with  moderate  hysteresis,  these  latter  effects 
can  be  treated  as  a  disturbance  and  targeted  with  a  feed¬ 
back  loop.  A  simple  control  design  utilizing  the  partial 
inverse  compensator  and  feedback  loop  is  depicted  in 
Figure  2.  We  note  that  this  design  was  used  to  establish 
the  potential  capabilities  of  the  partial  compensator  and 
is  not  optimal.  We  are  currently  investigating  the  place¬ 
ment  of  the  partial  compensator  inside  the  feedback  loop 
along  with  adaptive  temperature  compensation.  Exam¬ 
ples  illustrating  the  performance  of  the  partial  inverse 
compensator  are  provided  in  the  next  section. 


4.  Control  Implementation 

We  illustrate  in  this  section  the  implementation  of  the 
control  strategies  depicted  in  Figure  2  in  the  context  of  a 
high  speed  Terfenol-D  driven  milling  application  shown 
in  Figure  3.  This  represents  a  prototypical  application 
of  the  type  developed  at  Etrema  Products,  Inc.  The 
milling  specifications  include  a  3000  rpm  rotation  rate 
and  cutting  tolerances  of  approximately  ±1  micron. 

The  inverse  to  the  anhysteretic  curve  for  the  applica¬ 
tion  was  specified  to  be 


Ia(Xc) 


'  0.55x°-51  ,  0  <  xc  <  190 
<  0.0106xc  +  5.966  ,  190  <  xc  <  450 

0.45  x  10~4(xc  -  350)2-2  +  9.6  ,  xc  >  450 


to  accommodate  the  slight  degree  of  asymmetry  exhib¬ 
ited  by  the  data  (see  Figure  4).  This  represents  the 
inverse  of  the  general  anhysteretic  representation  in  (2) . 
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Figure  2.  Control  design  utilizing  the  partial  inverse  compensator  in  a  feedforward  loop.  A  feedback  loop  targets 
the  disturbance  due  to  hysteresis. 


The  performance  of  three  control  designs  is  summa¬ 
rized  in  Figure  4.  The  first  consists  of  an  open  loop 
control  law  specifying  the  commanded  position  with  no 
compensation  for  nonlinearities  or  hysteresis.  As  illus¬ 
trated  by  the  phase  plot  comparing  the  commanded  and 
measured  cutting  head  positions,  there  are  significant 
saturation  effects  and  hysteresis  which  then  produces 
errors  in  excess  of  70  microns  when  the  measured  po¬ 
sition  is  considered  as  a  function  of  time.  The  second 
design  utilizes  the  partial  inverse  compensator,  as  illus¬ 
trated  in  Figure  2,  but  without  the  feedback  loop.  From 
the  phase  plot  in  Figure  4b  it  is  observed  that  this  elim¬ 
inates  most  of  the  nonlinearity  due  to  saturation  but 
neglects  the  inherent  hysteresis.  Hence  there  is  a  slight 
phase  delay  in  the  time  history  of  the  measured  cutting 
head  position.  Figure  4c  illustrates  the  performance  of 
the  partial  inverse  compensator  when  the  feedback  loop 
is  included.  While  operating  under  a  slightly  different 
command  signal,  it  can  be  observed  that  the  feedback 
adequately  attenuates  the  disturbances  due  to  hysteresis 
for  the  operating  regime.  This  design  produced  a  cutting 
accuracy  of  ±15  microns  which  has  subsequently  been 
reduced  to  ±2  microns  using  a  narrowband  control. 

In  summary,  the  use  of  a  partial  inverse  compensator, 
based  on  the  anhysteretic  magnetization  or  polarization, 
provides  a  viable  design  tool  for  certain  hysteretic  smart 
material  systems  when  combined  with  a  feedback  loop 
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Figure  3.  High  speed  milling  application. 


which  targets  the  hysteresis.  The  advantage  of  this  tech¬ 
nique  is  the  efficiency  with  which  it  can  be  implemented 
and  the  observation  that  it  appears  less  sensitive  to  ac¬ 
curate  initial  conditions  than  the  full  inverse  compen¬ 
sator.  For  systems  which  exhibit  significant  hysteresis, 
however,  the  full  inverse  compensator  may  be  required  to 
attain  the  design  specifications  including  accuracy  and 
minimal  phase  delays. 
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